OUTP-01-21P 
April 2000 



On the Microscopic Spectra of the Massive 
Dirac Operator for Chiral Orthogonal and 
Chiral Symplectic Ensembles 



Frank Abild-Pedersen (1 ) and Graziano Vernizzi^ 2 ^ 

(1) The Niels Bohr Institute 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
e-mail: horsboll@alf.nbi.dk 
and 

(2) Department of Theoretical Physics, Oxford University 
1 Keble Road, Oxford, 0X1 3NP, United Kingdom 
e-mail: vernizzi@thphys. ox. ac.uk 



Abstract 

Using Random Matrix Theory we set out to compute the microscopic correlators of the 
Euclidean Dirac operator in four dimensions. In particular we consider: the chiral Or- 
thogonal Ensemble (%OE), corresponding to a Yang-Mills theory with two colors and 
fermions in the fundamental representation, and the chiral Symplectic Ensemble (%SE), 
corresponding to any number of colors and fermions in the adjoint representation. In both 
cases we deal with an arbitrary number of massive fermions. We use a recent method 
proposed by H. Widom for deriving closed formulas for the scalar kernels from which all 
spectral correlation functions of the %GOE and %GSE can be determined. Moreover, we 
obtain complete analytic expressions of such correlators in the double microscopic limit, 
extending previously known results of four-dimensional QCD at (3 = 1 and (3 = 4 to the 
general case with Nf flavors, with arbitrary quark masses and arbitrary topological charge. 



1 Introduction 



One of the most successful and well-established physical applications of Random Matrix 
Theory (RMT) is the analysis of Quantum Chromodynamics (QCD) at low energies. In 
particular, the spectral statistical properties of the Euclidean Dirac operator in the infrared 
regime, can effectively be described within a RMT approach. In fact, there is a close 
correspondence between the finite-volume partition function of four dimensional QCD in 
the low-energy limit and the partition function of a RMT with the same global symmetries 
[0 [@] 0- This theory is called chiral RMT (%RMT), because of its chiral content. The 
four-dimensional QCD Dirac operator with Nf fundamental fermions and gauge group 
SU(N C = 2) or SU (N c > 2), is described by the chiral Orthogonal (xOE) or chiral Unitary 
ensemble (xUE) , respectively. The case of adjoint fermions and SU(N C > 2) gauge group, 
corresponds to the chiral Symplectic ensemble (xSE) [||]. These agreements are valid only 
in the microscopic limit of ;\TiMT, in which the universality-classes of the three chiral 
random matrix ensembles, manifest explicitly ||. In this particular limit, one "magnifies" 
the Dirac operator spectra around the origin (zero virtuality) which, through the Banks- 
Casher relation ||, reflects the existence of a non- vanishing expectation value of the chiral 
condensate £ = (ipip), i.e. a chiral symmetry breaking at low energies. Indeed, Leutwyler 
and Smilga showed |7| that when S ^ suitable spectral sum-rules must hold. Such 
sum-rules can be written as integrals over the microscopic spectral density of the Dirac 
operator given by 

P.{*)= lim t^;P(t^;), (1-1) 

where V4 is the Euclidean space-time volume. The universality of %RMT in the micro- 
scopic limit (and therefore the universality of p s {z)) supports the idea that the QCD Dirac 
spectra are universal in the large- volume scaling limit [|TJ || ||. Several facts confirms 
this scenario, such as analytical results and extensive universality studies of microscopic 
spectral correlators JlO| [llj |12 Jl3| |14| as well as good agreement with numerical 



results obtained in Lattice QCD simulations for massless |D| and massive \Tfi\ fermions. 



For more details and a complete list of references we suggest some excellent reviews [|17 



in which it is also possible to find information and references about non-chiral RMT 
and the corresponding three-dimensional QCD theory. 

In the fundamental paper |J the large- N behaviour of the orthogonal polynomials relevant 
for the chiral Unitary ensemble is proven to be universal near the origin, in the scaling 
limit x = N 2 X, where x is kept fixed. From this it follows that all microscopic correlators 
have the same universal behaviour. Also in the so-called double microscopic limit of ;\UE 
where both masses and eigenvalues are scaled at large- N, one sees a similar universal be- 



1 



haviour of the orthogonal polynomials associated with Nj massive flavors. Consequently 

In the xOE and xSE the question 
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all microscopic correlators are seen to be universal 
whether the microscopic correlators in general are universal or not is still open. By relating 
the massless kernels of xOE and xSE to the xRMT with complex elements ((3 = 2), which 
is known to be universal, it is shown in ref. ||13|| ||14j| that microscopic universality persists 
under certain smoothness assumptions. Though a similar universal behaviour in the dou- 
ble microscopic limit for (3=1 and (3 = 4 remains to be proven, it is widely recognized 
that universality in these two models is a reliable consequence as well. 
So far several results have been explicitly obtained for the case of massless flavors. For 
instance, the microscopic spectral density is known for all chiral ensembles in the case 
of an arbitrary number of massless flavors [|J || |19| [20|. Also, the cases concerning an 



arbitrary number of massless flavors and an even number of massive flavors in non-xUE 
have been derived in [pl| and [^1|] respectively, and recently the cases of non-xOE and 



non-^SE have been solved in 

In the case of massive flavors, the situation is as following. 

m 



spectral correlators of the xUE is given in [T 



A general solution for the 
However, computing the microscopic 



spectral density in the massive (non-)xOE and massive (non-)xSE poses many difficulties 
arising from the use of skew-orthogonal polynomials and their behaviour at large- N. In 
p4[ p5| these ensembles have been studied and a solution is presented for doubly and (3- 



fold degenerate masses. Also, the general result of |12| could in principle give the spectral 
density for all the three ensembles with an arbitrary number of flavors and masses, but it 
seems technically difficult. 

Recently, H. Widom introduced a new technique in [2|| for dealing with OE and SE us- 
ing the standard orthogonal polynomials of the UE and thereby avoiding skew-orthogonal 
polynomials. This technique has already been used successfully in [22[] for studying the 
massless non-xOE and non-^SE. In this paper we apply the same technique to the %OE 
and xSE with an arbirtary number of massive flavors. Applying this method presents sev- 
eral advantages. First of all this approach seems to be a promising alternative to the usual 
one with skew-orthogonal polynomials, and it gives a different point of view of the general 
problem of determining the scalar kernels of the OE and SE. Secondly, since this technique 
deals with standard orthogonal polynomials only, then it is possible in principle to get the 
microscopic limit of all the correlation functions from the already well-established results 
on the universal microscopic limit of orthogonal polynomials in the UE. Finally, the or- 
thogonal polynomials for the general massive case naturally lead to a unifying notation 
which is very helpful when used with this new technique. 

This paper is organized as follows. In the next Section we shortly address the definition 
of the chiral matrix model which is relevant for four-dimensional QCD. In Section 3 we 
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show how it is possible to compute correlation functions for the x^E and xOE by means 
of the technique of Widom. Such a technique requires as a basic ingredient, the explicit 
evaluation of orthonorma/ polynomials for the general massive case. These are explicitly 
determined in Section 4 (with their normalization factors) in a closed formQ. Within this 
formalism, the degenerate massive, and the massless case are obtained as particular cases 
of our general formulas. Moreover, we derive some equalities among orthogonal polyno- 
mials with different number of flavors. In Section 5 we apply the machinery described in 
Section 3 to the case of interest here, and finally in Section 6 we derive an expression for 
the scalar kernels of the massive %OE and massive %SE. The microscopic limit of our final 
formulas is discussed in Section 7. 

Let us finally remark that in this paper we do not address the question of universality 
(which is actually considered here as a working hypothesis) of our results as well as ex- 
plicit numerical calculations. We postpone both of them, and in the following we focus on 
the application of the new technique by Widom to the massive case of four-dimensional 
QCD. 



2 Chiral ensembles 

A chiral random matrix model with the same symmetries of the QCD partition function, 
can be set up by replacing the Dirac operator V with a suitable constant off-diagonal 
block random matrix. The off-diagonal structure stems from the anticommutation relation 
between D and 75. Thus xRMT for QCD in four dimensions is defined by the partition 
function U [§, 

N f 

zW\{m f }) = JdW\{ det (V + m f ) e - N ^ v ^w) ? (21) 




» = ,w. „ • (2-2) 



with 



where the Dyson index (3 = 1,2,4 labels the different chiral ensembles. The rectangular 
matrix W of size Nx(N+\v\) is real, complex or quaternion real for (3 = 1,2, 4, respectively. 
The integration measure dW is the Haar measure of the group under consideration, the 
integer v is related to the topological charge @ and 2N + \u\ is the space-time volume. 
The polynomial potential V in the exponential can be replaced by a pure quadratic term 

1 The same set of polynomials have already been obtained in llOl, but only in an recursive form. 
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V(WtW) = W^W because of the universality of the matrix model in the microscopic 
large- N limit |J jl3R . Throughout this paper we will deal with an arbitrary topological 
charge v. 

The matrix integral can be rewritten, up to an irrelevant overall constant factor, in terms 
of the eigenvalues \ of the Hermitian positive-definite matrix W^W, that is: 

N ( N f \ 

ZW\{m f }) = fn ^n^ + ^Af^V^ |A({A fc })| 

N 

n^A i «;(A i )|A({A fe })| /3 , (2.3) 

i=l 

where A({Afc}) = Ilm>n(Am — A n ) is the Vandermonde determinant and 

w (\) = J] (A + m^X^+^-'e-^ (2.4) 
/=i 

is the weight function on I = [0, oo[. These two functions define all the spectral statistical 
properties of the matrix model in eq. ( |2.3| ). The corresponding properties for the model 
in eq. (|2.1| ) are expressed in terms of the real eigenvalues & of the Dirac operator T>, by 
means of = A j. 

3 Spectral correlation functions for (3=1 and (3 = 4 

The m-point spectral correlation function for all three chiral ensembles is 

/V' m 

^(Ax,..,^) = J (YlTrSfr-WiW)) (3.1) 
= det [IffCAi.A,-)] , (3.2) 

l<t,j<m 

where the expectation value in the first line, is understood with respect to the partition 
function in eq. (|2.3j ) and the second line follows from a well-known result |27] which 



expresses all the correlation functions in terms of the determinant of the kernel of suitable 
polynomials^]. For (3 = 2 such polynomials are orthogonal polynomials w.r.t. the weight 

2 For [3 = 1 and /? = 4 the determinant in eq. ( |3.2| ) is understood as the quaternion determinant of a 
matrix kernel (for details see p7[). 
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w(x) in eq. ( j2.4| ), whereas when (3 = 1 and (3 = 4 one usually chooses skew-orthogonal 
polynomials w.r.t. w(x). Indeed, one does not need to introduce skew-orthonormality 
since the matrix kernel is largely independent of the particular choice of the polynomials 



f29| . Therefore, if one defines 



<p 3 (x) = Pj (x)^w(x), j = 0,l,...,2iV-l, for/? = 4, (3.3) 
ifjix) = Pj {x)w{x), j=0,l,...,iV-l, for/3 = l, (3.4) 

where Pj{x) are arbitrary polynomials of order j, then the matrix kernels can be written 
as pOfl : 

iWf~ „,\ c( 4 ) 



„\ - ( s n(?,v) S n D(x,y) 
IStffay) S$\y,x) 



and 



K$(x,y)=( m ^f(^) ) , (3.6) 

ISy(x } y) -e{x-y) S ( N '(y,x) 

where e(x) = sgn(x)/2, S^\x,y) equals the scalar kernel of a suitable operator Sjy , and 
ISx\x,y) and DS^\x,y) are the kernels of the operators Tsffi, DSffl with 7, D identi- 
fied as the integration and differentiation operators, respectively. All the matrix elements 
in eq. ( |3.5[ ) and eq. ( |3.6| ) are determined once the scalar kernel S^\x,y) is known. 
Indeed, the scalar kernel S^\x,y) can be expressed in terms of ipj functions only and it 
is a matter of fact that the choice of skew-orthogonal polynomials leads to the simplest 
possible expressions (see for instance |2"7fl). However, skew-orthogonal polynomials are 



difficult to determine and furthermore they seem unsuitable for straightforward manipu- 
lations at large- N as opposed to the standard orthogonal polynomials^. This fact has led 
some authors to use a different approach |p6]| . This approach uses orthogonal polynomials 



only and, without introducing skew-orthogonal polynomials, H. Widom derived general 
formulas for the kernels S^\x,y) expressing it as corrections in addition to the unitary 
kernel K^\x,y), with N replaced by 2N when (3 = 4 and w replaced by w 2 when (3 = 1. 
Because all quantities are expressed in terms of orthogonal polynomials, calculations and 
large- N asymptotics are straightforward in this approach. 

In this paper we adopt the method developed by Widom in order to determine the func- 
tions Sjpfe, y) and through that we analyze the original physical problem. In p2| the same 



technique has already been applied successfully to the case of three-dimensional QCD with 



3 A formula useful for studying the asymptotics of skew-orthogonal polynomials, has been derived 
recently in [RlJ. 
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Nf massless fermions. In the remaining part of this section we briefly sketch this technique. 
For more details we refer the reader to the original works [26] [29] and to \2% . 



The first step of this technique, is to build-up a Hilbert space Tt from the functions ipj 
defined as: 

(3.7) 



0,1, 



N-l 



iPj[X) =Pj{X) yjW\X) , J 

According to eq. fl3.3| ) or eq. ( |3.4[ ), in what follows we always assume that N should be 
replaced by 2N when j3 — 4 and w should be replaced by w 2 when (3 = 1 (if not explicitly 
stated otherwise). Under suitable hypotheses on the measure^ w(x) = e~ v{ ~ x \ TC is defined 
as the linear space spanned by the functions (fo, yj\r-i- In this picture the scalar kernel 
Kffl(x, y) can be considered as the kernel of the projection operator^ K onto the Hilbert 
space H, and it may be written as 



K$(x,y) = <Pj(x)<Pj(y) 



N-l 



j=0 



x-y 



1 

1 



(3.f 



where = kN-i/kjy, with k^ identified as the highest coefficient in pj^(x). The kernel 
Sx\x,y) can be written as the unitary scalar kernel K^\x,y) plus extra terms, whose 
number is independent of N and closely related to the commutator [D, K\. The kernel of 



the commutator [D, K] is |26 
[D,K]{x,y) = 



dx dy 
aN (<Pn, <Pn-i) x 



C(x)-C(y) A{x)-A(y) 
(y) 



A(x X )-A(y) 



B(xT- 



x-y 



x-y 




where 

A(x) = -A N (x) 



-V'(x), B(x) = B N (x), and C(x) = -^B N ^(x) 
2 cln-i 



are rational functions with the same poles as the ratio w'(x)/w(x), and 

r+oo 

A N (x) = a N if N {z)if N _ 1 (z)U{x,z)dz, 
Jo 

r+°° 

B N (x) = a N / ^>n{z) U(x, z) dz , 
Jo 



(3.9) 

(3.10) 

(3.11) 
(3.12) 



4 In particular w'{x)/w{x) must be a rational function over I. 

5 For that it is crucial that the polynomials {pj(x)} are orthonormaZ. 
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with 



and 



V'(r) — V'(z) 
U(x, z) = { > , V'(x) = X (P)V'(x) , 



x — z 



1 for/? = 4 

2 for = 1 ' 



(3.13) 



(3.14) 



stemming from the change in notation. Let and n Xi denote the pole orders of w'(x)/w(x) 
at infinity and Xi, respectively. Then the functions 



X (f N -l 
(x -Xj) _fe_ ViV-l 



x k (p N (0 < k < Uoo) , 
(x-Xi)- k - x ip N (0 < k < n x .) , 



(3.15) 
(3.16) 



span a subspace 7"4ub C and a corresponding subspace C Ti -1 both of dimension 
n = "^oo + J2i n xi- The subspace 7Y su b is determined by the condition that 7"4ub C 7Y, 
whereas the subspace 7"^ b is determined by the n orthogonality conditions set up by the 
functions 

(f N - k (k < rioo) ; ip k (kKn-n^) . (3.17) 

Let ipi, ip n G Tt and ip n +i, ■■■,i J 2n £ / ^" L denote these 2n linearly independent functions, 
then one has 



2n 



[D,K](x,y) = ]T 1 p i (x)A ij t/; j (y) , 



(3.18) 



and once the ^'s are chosen, the symmetric constant matrix A = [Aij] is uniquely deter- 
mined. Furthermore, the matrix A is always in a block off-diagonal form: 



A^ = i,j<n or i,j>n 



One also defines the matrices: 



jf jf dxcfy e(x-y)i/ji (x) i/jj (y) , 
\ 



(3.19) 



(3.20) 



J 

C 



\ nxn 



2nx2n 



(3.21) 
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with O nxn , Inxn being the n x n null matrix and the identity matrix, respectively. The 
operator e is defined as 

(eif>i)(x) = [ dxe(x - y)tpi(x) . (3.22) 



i 



Defining A , C and C 00 as the matrices obtained by deleting from the corresponding 
matrices the last n columns, the last n rows and the last n rows and columns respectively, 



one finally has the main result of [26|, that is 

2n 

S$\x,y) = K%\x,y)- £ [AC (I - BACy^^e^iy) , (3.23) 

i<n ,j=l 
2n 

S$\x,y) = K%\x,y)- £ [A^C^^e^y) , (3.24) 

i>n,j=l 



remembering once again that iV — > 2N in eq. ( |3.24|) and w — > w 2 in eq. ( 3.23| ). 



A remarkable feature of eq. (|3.23j ) and eq. (|3.24|) is that the scalar kernels S^\x,y) 
naturally are expressed as a finite number of corrections to the unitary kernel K N (x, y), 
a number which is independent of N. This fact will prove its importance in Section |7] 
where we study the microscopic limit of the quantities of interest here. Finally, inserting 
eq. ( p.23| ) and eq. ( |3.24| ) into eq. ( |3.6| ) and eq. (|3.5|), one can calculate eq. ( |3.2| ) and 



through that all possible spectral correlation functions defined by eq. (3.1) 



4 Orthogonal polynomials for Nj massive fermions 

In this section we apply the technique of Widom, described above, to the general mas- 
sive chiral case with Nf flavors. The first ingredient we are looking for are polynomials 
Pn f (x; mi, • • • ^ m N f ) orthonormal with respect to the weight function 

N f 

w {N f' a \x) = l[(x + m\) x a e~ x (4.1) 
i=i 

defined on the interval X = [0,+oo[, where a is a real non-negative constant. In this 
formalism, it is understood that w^ ' a \x) = x a e~ x . The weight in eq. ( (4.1|) is seen to 
be slightly different from the one in eq. (|2.4| ), but if we identify a = j3{v + l)/2 — 1 and 
c = N/3/2, then we observe that 
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the variable a could be negative (e.g. a = — 1/2 for (3 = 1 and u = 0). Although we 
suppose a > throughout this paper, the final results are valid also for a > — 1 i>za 
analytic continuation PBfl; 



• the exponential function in eq. fl4.1p should be e~ cx indeed. Here we suppose c = 1 
instead, but at the end of Section || we shall extend all the results to the general case 
c 1, by means of scaling arguments. 

(Nt a) 

In order to determine the polynomials P n explicitly, we fix some useful notation at 
first. We define the scalar product of two real functions f(x), g(x) with respect to the 
weight function w( N f' a '(x) by 



(f,g) Nf>a = j^dxw^ a \x)f(x)g(x) , (4.2) 

and furthermore we shall use also the following shortened and suggestive notation: 
Pi N f' a \x;m) = P^ a \x; mi ,...,m Nf 



Pi N f> a \x;m 1 ,...,m Nf - 1 ,0) , 
Pf^^ijmi, . . . ,mj_i,m i+ i, . . .,m Nf ) . 



Moreover, let k^ f '°^ denote the highest coefficient in P^ f '°^ 



kn x n + . . ., and let us 



define the (N f + 1) x (N f + 1) matrix 

/ P(°> a) (x) 



P^(-ml) 

p(0,a) 



-m 



p(0,a) ( o ^ 
r n+l I ,n N f , 



r n+Nj \ m 



P 



(0,a) 



n+Nf 1 



-m 



(4.3) 



Finally, we identify A^/' a \ i = 1, . . . , Nf + 1, as the Nf x Nf sub-matrix obtained by 
omitting the first row and the z-th column from An f ' a \x) in the definition fl4.3|) . Notice 
that all A matrices are expressed in terms of orthonormal polynomials P^°' a \x) with 
Nf = flavors only. Since the polynomials P^ ,a \x) are orthonormal w.r.t. the measure 
x), then they are necessarily proportional to the generalized Laguerre polynomials 



w 



(0,o), 



LW(x) @,i.e. 



if)(x)=LW(x)/^, ^= r(n + ? +1) 
v n! 



(4.4) 
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Therefore, from L^\x) = (— \) n x n /n\ + . . . one reads the highest coefficients 

fr(0,«) _ 

r<J^ — 



(4.5) 

T(n + a + 1) n\ 
Now, it is well-known that the orthonormality condition 

(Pj Nf ' a \x; m), pj Nf ' a \x; m)) Nf , a = h 3 (4.6) 

uniquely determines polynomials P^ Nf ' a \ up to a relative sign. Under the hypothesis that 
all the masses {m,} are distinct, these polynomials can, according to Christoffel's theorem 
3], be represented in terms of P^' a '(x) as follows: 



p( N f> a \x;m) 



1 det 



mi) 



(4.7) 



The normalization factor h^n f,a ^ is 



( m ) 

and the highest coefficient is 



■l)^^Tdet 

fvn 



.(N f ,a) .(JV>,a) ' 



(4- 



fc (iV/,a) ( 



111 ) 



\ 



;-i)^a4 m A; (0,Q!) 



det 


a(^/.«) 1 
iv n,AT / +l 


det 


r A (^/,")i 



(4.9) 



where fc^ ' a ) is the coefficient given in eq. (|4.5|) . For an explicit derivation of eq. 
and (f4.9|) we refer to appendix Eq. ( |4.7| ) is completely symmetric under permutations 
of the masses: although the determinant in the numerator is completely antisymmet- 
ric, the algebraic square root in the denominator is also completely antisymmetric (with 
— 1). This observation is consistent with the fact that also the weight function 



-1 



in eq. (4.1) is completely symmetric under permutations of the masses. Moreover, the 
case of a degenerate mass is understood in a limit sense of eq. (|4.7Q . For instance, 
if there are only two degenerate masses m\ = m 2 , we just take the limit of eq. (|4.7| ) as 
rrii. Dividing both numerator and denominator by m\ — mf, and substituting the 



m 2 



third row {P}° (— mj)} with the combination {Pi U ' a> (— m%) — P^ u ' a> (— w^)}, in all the 



(0,a), 



)(0,a). 
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determinants, we obtain in the limit m\ — > rri2, that the third row effectively is replaced 
by the first-order derivatives dP^ ,a \y) / dy at y = —m\. In the case of a degenerate mass 
rrik of multiplicity I, I > 1, we replace in all the determinants of eq. ( j4.7| ) the corresponding 
rows k + 1, k + 2, . . . , k + 1 with the derivatives of order 0, 1, 2, — 1 of the polynomials 
Pf' a \y) evaluated at y = —m\. It is worthwhile to remind here that degenerate masses 



are definitely needed for the (3 = 1 case, because w — > w 2 effectively implies Nf 
i.e. {m 1 ,m 2 , . . .,m Nf } -> {m u m u m 2 , m 2 , . . .,m Nf ,m Nf }. 



2N 



Let us give two explicit examples where eq. ( [Of ) is applied. In the case Nf 
reads 



1, it 



n\(n + 1)! 



\ + a + l)r(n + a + 2) 



(x + m\)\Jhn ,a \mi) 



where 



e a) K) 



n! 



T(n + a + 2 
In the case Afy = 2, with degenerate masses m\ 



LtU-rnl)L^{-ml) 
m 2 = m, it reads 



(4.10) 
(4.11) 



det 



L 



(a) 
n+2 



X 



\ 



p(2,a). 



x; m, m) 



4 a) (-™ 2 ) 4 a i(-™ 2 ) 4 a i(-™ 2 ) 

Ui a) '(-^ 2 ) 4ti(-™ 2 ) L%{-m?)) 



(x + m 



2\2 



IL n IL n+l' i n+2 ,i n 



(2.«)/ 



(4.12) 



and 



n!(n + 1)! det 



m, m 



4 Q i(-™ 2 ) 

42i(-™ 2 ) 




LW(-m 2 ) 
L^'i-m 2 ) 



L ( n a U-m 2 ) 



T(n + a + 2)T(n + a + 3) 



(4-13) 

All the derivatives can easily be evaluated using the property dL^{x)/dx = — L^-l, n > 
(iteratively, if needed). 

These examples are just two particular cases of a general fact: all the orthogonal polyno- 
mials for the general massive case are expressed nicely in terms of generalized Laguerre 
polynomials by means of eq. (|4.7| ) , ( |4.8|) and (|4X 



Now, for the sake of future use we point out two remarkable properties of the polynomials 



in 



(N f ,a) 
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1. Let us consider the set {mi, . . . , , 0} consisting of Nf + 1 flavors. In this case eq. 
Q4.6Q reads 

% = (^ /+1,a) (x;m,0),Pf' +1 ' a) (x;m,0)> JV/+lia (4.14) 

■> 1 m Nf+1 =0 

that is, the polynomials Pn f+1 ' a \x; m, 0) are orthonormal with respect to the mea- 
sure x a+1 ]li=fi (x+mf)e~ x dx. On the other hand, also the polynomials pi Nf ' a+1 ^ [x; m) 
are orthonormal w.r.t the same measure, i.e. 

= (Pf />a+1) (x; m), pf m)) W/ , Q+1 . (4.15) 
From uniqueness of orthogonal polynomials w.r.t. a given measure it follows 

p(N t +lv)( x . m 0) = p(N f 0+l)( x . m ) ) (4.16) 

for every n and iVy|] 



2. Christoffel's theorem eq. ( |4.7|) can also be stated as 

det 



P!- N f> a Hx: m) = d 



(N f ,a) t 



r n 



(N f -l,a) 



r u+l {X, "l^j 



III I 



-m; 



2. 



fx + m. 2 



for i 



(4-17) 

iV/, where the coefficient c^/'^ is easily determined by comparing the 



highest coefficients on the two sides of eq. ( |4. 17] ) , i.e. 

k { n N i' a \m) 1 



SN f ,a) 



111 ) 



^r 1,a V^)^" 1,a) (-^ 2 ;m^ 



(4.18) 



And, by means of the Christoffel-Darboux formula |33[ , one can equivalently write 
eq. (PTTD as 



(JV/.a), 



m 



J=0 



(iV>-l,a) 



-m^m^P, 



(N f -l,a) 



x;m^) , (4.19) 



6 A trivial check of formula ( 4.1 6| ) is for Nf = 0. Indeed, if we put mi = in eq. (4.10) then 
we obtain exactly Pn°' a+1 \x), eq. (4.4). Similarly, if we put m = in eq. pl^) then we obtain 



P^foO.O) = p£' a+1) 0) = P^ a ^(x) 
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for i — 1, . . . , Nf where the coefficient d^/ ' a \m) is determined as before by compar- 
ing the highest coefficients in eq. ( |4.19| ), i.e. 



JN f ,a) 



dn,i ) ( m ) — ,^N f -i,a), > n(ff/-M)/ 2. - ' ( -1-20) 



The coefficients fen^'^'s which appear both in eq. ( 4.18|) and in eq. (|4.20|) , are given 
explicitly in eq. (|4.9|). 



These two considerations will be useful in the evaluation of some integrals appearing in 
next Section. 



5 The Hilbert space 7i 

Once the orthonormal polynomials are known then the next step is to build-up a suitable 
Hilbert space according to the technique described in Section |3|. In order to do that, from 
the orthonormal polynomials Pn f (x; m) we define the functions (Pj{x) as in eq. ( |3.7| ): 

<pj(x) = Pj Nf,a) (x; m)^Jw( N f^(x) , j = 0, 1, . . . . (5.1) 

The Hilbert space is Ti = span{y9 , <fi, ■ ■ ■ , <Pn-i}, and the kernel of the projection operator 
K onto Ti and the kernel of the operator [D, K] are given in eq. ( |3.8|) and eq. (|3.9|) , 
respectively. In particular, we notice that the former can be nicely written in a very 
compact form: 

N-l 



K%\x,y) = Yl Pj N " a \x; m)P^ Nf ' a \y; m)^/w( N f^{x)w( N f^(y) (5.2) 



3=0 

a N-l,N f +l\ m ' V V) 



yJ w (Nf,<*)(x)w( N f< a ){y) , (5.3) 



where we used eq. ( |4.19[ ) in the second equality. 
Since w^ N f' a \x) = i.e. 



N -l 
»=i 



V(x) = x — a log x — \og(x + ml) , (5.4) 
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eq. ( |3.13| ) reads 



U{x,z)= X {/3) 



xz i^i ( x + m i ) ( z + m i ) 



(5.5) 



where is defined in eq. ( |3.14| ). Therefore, substituting eq. ( |5.5| ) and ( p.l| ) into eq. 
( |3.11|) and eq. ( [3.12| ) we obtain the matrix elements in eq. ( |3.9|) : 



(5.6) 

(5-7) 
(5.8) 



A(x) 

X 


-y 


= a N x(f3) 


xy V 


B(x) 


-B(y) 


B N (x) - 


-B N (y) 


X 


- y 


x - 


- y 


C{x) 


-C(y) 


a N B N _i(x) — B N 


X 


- y 


a N-l 


x-y 



-a N x((3) 



Nf AN,N-1) i_ 

1 \ , __2 aajy 

2a N J t^i( x + m i)(y + m i) 

j(N,N) 



oI ( N > N ) N f 

+ 



xy 



{ (x + mf){y + mf_ 



where a at = k^£*\m) /k^ f,ull {ja\) and the constants I . /, 



S N f> a )i 



' A '' ' Ii N '^ are defined as 



/P ( N f> a ) f_i \ 



(Pi 



(N f ,a) 



Z 



N 



z + mf 



/N f ,a 



i = l,2,...,Nf 



(5.9) 
(5.10) 



The above integrals are evaluated in appendix [B] using the previously stated properties of 
orthonormal polynomials for general Nf. Thus we have 



T (N,N-1) 
T (N,N) 
T (N,N-1) 



4^;ir 1) (m,o)pf_r ) (o;m) 

(N f ,a) 



jN f +l,a~l) 
a N,N f +l 



d 



(N f ,a) 



N.i 



m)P, 



(m,0)P- / '" ; (0;m) 

(N f ,a) 



N-l 



T (N,N) 

which appear in eq. ( |5.6|) , (|5.7|) and ( |5.8| ). 



—m^ m) 



—m i ; rri j 



Now we can construct the ipi(x) functions. Since the logarithmic derivative of the weight 
is a rational function 



d w (Nf><*)( x )/dx _ a + ^ 



w 



( N f> a )(x) X 



i=0 



x + mi 



- 1 



(5-11) 
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with no poles at infinity and simple poles at x = (for a ^ 0) and at x = —mf, i = 
1, . . . ,Nf, we have that the subspace Ti su b has dimension n = Nf + 1. The number of 
linearly independent functions ipk 6 Ti we are looking for is thus equal to n = Nf + 1. 
According to eq. ( |3.15|) and ( 3.1 6|) we choose the following linear combinations: 

in {N f' a) (r) in {Nf ' a) (r) 

if N {X <f>N~l \ X ) -i AT , 1 (C 10\ 

ipi{x) = ■ — 2— + A ; — o - ' « = 1, + 1 , (5.12) 

where we used the shortened notation Oi = rrii for i = 1, . . . , Nf and 0jv>+i = 0. We 
determine all the coefficients Dj by requiring that the functions ipi are non singular at 
x = —erf, since they have to belong to the Hilbert space Ti. This means that the numerator 
in eq. ( |5.12j ) should be vanishing at x — —erf, i.e. 

Di = - " } *' } , i=l,...N f + l . (5.13) 

Thus all the functions ipi G Ti are determined completely. In appendix we obtain the 
following compact final expression: 

(N f +l,a) 



Mx) = fc^^M ^- 1 ; ;m '^V wW,a) W , i = l,...,JV/ + l. (5.14) 

fckr -i ' (m,(7i) 



We now look for the functions ipN f +2, ■ ■ ■ ,4>2N f +2 £ Ti 1 - . We choose again to consider 
the combination 

V>jv /+ i +i (x) = ■ — — + D i ■ —, I = 1, . . .,Nf + 1 , (5.15) 

J x + of x + of 

where D/- are coefficients to be determined. As discussed in Section | the space S = 
span{</>o, • • • ,fN f } is a n-dimensional subspace of Ti, henceQ we can fix all the D^s by 
imposing that each ipN f +i+i G Ti is orthogonal to all the {<^ , . . . , <fN f }, i-e. 



ip Nf+1+i {x) <^j(x) = , (5.16) 

for i — 1, . . . , Nf + 1 and j = 0, . . . , iVy. When z = 1, . . . , Nf, by substituting eq. (|5.15| ) 
and ( |5.1| ) in the constraint equation ( 5.16p , we see that all the integrals are of the form eq. 

7 This fact implies N > n — Nf + 1, because Ti is TV-dimensional. 
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( |5.9|) or ( |5.10| ), which have already been solved explicitly in eq. ( |B.4| ) and ( |B.5| ). Therefore, 
the constraint conditions (|5.16|) can be written as: 

n _L a N,i \ In ) ■ 1 AT /r 17 \ 

A =- j Nf , a) . - . i = l,...,JV), , (5.17) 

( m ) 

with iV > iVj + 1. Notice that the bound N > iVy + 1 is fully consistent with our purpose 
of taking the microscopic large- N limit (see Section [?p. When i = Nf + 1 we have 



D ± _ <W'i~ (m,0) 



,(JV/+l,a-l) 
l N,Nf+l 

In appendix |C] we obtain a more compact and explicit expression for ^j^+i+i, i-e. 



WN f +i+i(x) - k N (mj {N l a) j— — 27 , (5.19j 



for 2 = 1 , . . . , Nf and 



.( K ,, a ),- Pf^ 1) (x;m) v ^^ 

1p2N s +2{X) - % l m J ^JV/.a-l)/ ; ~ ' (D.AIJ 



for z = JV/ + 1. 



6 The kernels (x,y) 



Now we have all the necessary ingredients to determine the symmetric 2n x 2n matrix A, 
which is defined by eq. (|3.18 ). In [26 it is proven that the matrix A always has the block 
form 

(O ■■ A \ 



A 



\ 4xn : U nxn J 2n 2n 

where nx „ is the n x n null matrix. With our choice of the ip functions, the n x n matrix 
A is always diagonal. Indeed, it is sufficient to prove that the off-diagonal terms A\j are 



(6.1) 
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identically zero. In fact, if such terms were not zero then they would be responsible for 
having mixed terms in [D,K](x,y), eq. fl3.18|) , of the type ^i(x)ipj^i(y) , that is mixed 
terms of the type l/[x(y + m?)] or l/[(x + mj){y + m?)] (i ^ j). But such mixed terms 
are definitely not there in eq. (|5.6|) , (|5.7Q and ( |5.8|) and therefore, through eq. ( |3.9|) , they 
do not appear at all in [D, K](x,y). Finally, from eq. ( |3.18 ), we necessarily conclude that 
A i:j = for i ^J. 

Let us define A4 = An as the diagonal elements of A, i = 1, . . . , Nf + 1. These elements 
are uniquely determined by comparing eq. (|3.9|) and eq. (|3.18|) : indeed, in our case it is 



sufficient to compare the highest coefficients in x and y. In particular, in eq. (|3.9| ) the 
terms proportional to P^ fi ' a \x)P^ f ' a \y) are contained in 



a N P^ a \x) C{x) - C{y) P^\y) 



x-y 



(6.2) 



whereas in eq. ( |3.18| ) and ( |5.12j ) such terms are contained in 



2P 



(Nf, a) 



N 



X) 



E 



A, 



1 (x + m|)(y + mf) 



+ 



A 



xy 



P 



(N f ,a) 



N 



(y) 



(the factor 2 comes from the symmetricity of A). Comparing eq. ( |6.3|) with eq. (|6l 
( |5.8| ) term by term we identify 



2A, 



(JV-l.JV-1) 



(Nf,a) 
N-1A 1 



m)P 



(N f ,a). 



AT-1 



-mf; m) 



2A 



Nf+l 



2 , n\ T (N-l,N-l) 



(6.3) 



and 



.,N f 
(6.4) 



from which, the matrix A is completely determined^. Moreover, it is a result of Widom 
that the matrix A is uniquely determined by the given choice of ipi functions [2f| . Therefore 
it is garantueed that all the lower-order terms in eq. ( |3.9[ ) and eq. ( |3.18| ) match exactly. 
To calculate the corrections, we need the antisymmetric matrix B in eq. ( |3.20|) . Unfortu- 
nately enough, the matrix elements (which are a total number of 2n(2n — l)/2 distinct 
elements) are much more difficult to evaluate than Ay, and therefore our final result for 



the corrections to the unitary kernel is given by eqs. ( |3.23p , ( |3.24| ), and ( |3.2(J| ), with -0's 
explicitly given in eqs. ( ^.14|) , ( ^.19| ) and ( |5.20|) and the matrix elements of A are given in 



Although such a final expression is not terribly simple, it is at least suitable for direct 



8 Notice that with our choice of ip functions, all the elements Ai are determined through the function 
C(x) only. 
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numerical evaluations. However, its major usefulness appears when one considers the mi- 
croscopic limit of it. In fact, one should remember that the physical interest is in the 
microscopic limit, since it is in this limit that universal properties of xRMT appear. In 
Section [7] we will study the kernel (x, y) exactly in this limit. 

Before doing that, let us mention that all the results of the last three sections can be 
generalized easily to the case where the weight function is 

N f 

w {N f> a ' c \x) = H(x + ml) x a e~ cx , (6.5) 

i=l 

i.e. where an additional parameter c > has been introduced in the exponential. We do 
not need here to repeat all the calculations above with this new weight function, because 
it is possible to write down the c-dependence for any quantity, explicitly, just using scaling 
arguments. For instance, from the identity 

N* 

poo 



1 = [°° dx [p( Nf,a \x; y/aa)] x a ]\{x + cm?) e~ x 
•' i=i 

N f 

poo . ,. r v _2 £ 

= c a+N f +1 dx[Pf f ' a \cx;V^m)]x a l[{x + m 2 i )e- cx , (6.6) 
•' Q i=i 

we can read off the orthonormal polynomials w.r.t. the new weight function in eq. ( |6.5| ): 

P- 1 '(x;m) =c 2 1 '{cx;\/cm) . (6.7) 

The same argument applies to all the other functions. In general, a quantity g c= i(x;m) 
defined with the measure «jW ,a ' (x)dx turns out to be related to the same quantity g c (x; m) 
defined with the new measure w^ Nf,a,c \x)dx according to this formula: 

g c (x; m) = c< ' g c=x [cx\ Van) , (6.8) 

where 7 is a suitable scaling exponent. The following table shows the scaling factor for the 
most important quantities in this paper. 
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Quantity at c = 1 


7 for c 7^ 1 


Quantity at c = 1 


7 for 1 




-iV/ - a 




1 


P\ Nf ' a \x-m) 


(a + N f + l)/2 




1 




n + (a + N f + l)/2 


C(x) 


1 


det[AT" a, (x)] 


(JV/ + l)(a+l)/2 




1 


^" a, (m) 


lV)(a + 2) 




3/2 




-(a + iV> + l)/2 


D D (L) 





43""(m) 


(l-a-JV»/2 


Aj 


-1 




1/2 


Bij 


1 


U(x,z) 


2 




1/2 


a N (m) 


-1 




1 



The only exception of this scaling rule, is given by the potential V(x) which, for c 7^ 1 is, 
V(x) = c — a hgx — J2kii( x + m D- 

We conclude this Section remembering once again that at the very end one has to perform 
an important substitution. Namely, it is a prescription of this method to replace N — > 2N 
when f3 — 4 and w — > w 2 for (3 = 1. The latter substitution leads one to consider double 
degenerate masses TOj, which means Nf — > 2iV/, and to take a — > 2a, c — > 2c. We 
summarize these rules here: 

/3 = 1 : iV/ -> 27V/, a -> 2a , c -> 2c , . 

/3 = 4 : iV^2iV ' {b ' 9} 

7 The microscopic limit 

The physical interest in studying RMT relies on its universal properties which appear in 
the large- N limit. In particular, we consider the double- microscopic limit of the scalar 
kernel: 

where Ci, C2, and /ij are kept fixed. In general, the microscopic limit of a quantity g c (x, m) 
in this paper is: 



g c (£,n) = lim A^<? c -i(^4,^) = ( ~V 1™ N^g^A—^ —fi) (7.2) 
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where the exponent 7 can be read off from the table at the end of Section 6 and 5 is a 
suitable exponent which is chosen such that a finite non-zero limit exists. According to 
eq. ( |7.2| ) we may obtain the microscopic limit of g c (x, m) simply by computing the mi- 
croscopic limit of g c= i(x, m) with the scaling x = (/N, m = /Lt/x/iV, and then evaluating 



the result at the points £ — > ( = ((3/2, fi — > fi = J/3/2 fi. Following this strategy when 
computing the microscopic limit of S^\x,y) in eq. ( |7. 1] ) , we consider c = 1 in the rest of 
this section, and then at the very end we substitute ( — > ( and /j, — ► /2. 
In eq. ( |3.23D and (|3.24 ) the scalar kernel S^\x, y) is expressed in terms of the unitary ker- 
nel K$ (x, y) and the ipi functions which all are determined explicitly in Section [5| Hence 
we can now compute the microscopic large- N limit of the scalar kernel (x, y). The basic 

(Nr a) 

ingredient is the set of the orthonormal polynomials P n and their microscopic limit. 
Such a limit is straightforward analytically. In fact, the polynomials are defined completely 
in terms of determinants of matrices of the form as in eq. (|4.3|) . If one naively takes the 
microscopic limit of such determinants, then it eventually ends up with an indeterminate 
form. The situation is actually analogous to the case of degenerate masses, i.e. one has to 
substitute rows or columns in the determinants with suitable linear combinations before 
taking the microscopic limit. For instance, the generic determinant can conveniently be 
written as (we use the shortened notation to = x,ti = —mf): 



det u [L^ +l (U)} det u [t[L^ +l) (U)] 



det[A JV / ' ] — det[P$$(ti)] - n v+.V; ,777 ~ n \+A , ,77 n a , 



with i, I = 0, . . . , Nf. In the last equality we iteratively used the recurrence relation 
xL^ +l \x) = (n + a + l)L^\x) - (n + l)L^li{x). Now substituting the microscopic 
scaling t,i — * Zi/N, one obtains the following large- A behaviour 

, N s N f {N f +l) (JV f + l)(c-iV f ) 

det[A^' a) ] ~ (-l)+f—N ' f a 1 det^ 2 J a+ i{2^)) , (7.3) 

il 

since limAr^oo L^\x/N)/N a = J Q (2y / x)/x Q / 2 . When the argument of the Bessel functions 
is complex (i.e., when ti = —mf), one will make use of J v {iz) = i u I v (z). Finally, from 
the previous simple arguments it follows that the microscopic limit of all the quantities 
evaluated in this paper, can be obtained using the determinant in eq. ( |7.3j ) instead of the 
determinant at finite- A given in eq. ( [4.3|) . 

Let us look at some explicit examples. For sake of simplicity, we consider Nf non- 

(Nt a) 

degenerate masses first. The microscopic limit of the normalization factor h N in 



eq. is: 

~ h (N f ,*)^ = [A(°V)] 2 , (7.4) 
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where fi = {fix, fi2, ■ ■ ■ , /J-N f } and 

A (a) (/ii,/i 2 , ...,HNf 



det[^- a I a+q (2fi 



p/i > 



(7.5) 



where p = 1, . . . , Nf, q = 0, . . . , iVf — 1. In this definition, we assume that the range of 
p and q is determined by the number of arguments of A^. We can also compute the 



microscopic limit of the highest coefficient k^'^'im) in an analogous way: 



( M ) = i , (7.6) 

because the ratio of the two determinants in eq. ( |4.9| ) equals unity in the microscopic large- 
N limit. From eq. (|4.7|) , ( |7.6| ) and ( |7.4|) we immediately read off the microscopic limit of 
the orthonormal polynomials for Nf massive fermions, that is: 

A< a) (£/x) 



p(iV/,a)/£. 



AH(M)nS(c + ^. : 



where (£, /it) is 



AS Q) (C,M)=det 



y ur a / a (2/ii) /il- Q / Q+ i(2/i 1 ) 



iv r — ex 

(-l)^C^^+iv / (2^) 



(7.7) 



The microscopic limit of the coefficients cf^f (m) in eq. ( fl.20|) is 



A (a) (^)n^j-M. 

(_l)*+iA(«)(/*) 
The microscopic limit of the weight w^ Nf ' a \x) is given by 



W 



(o = rn(c+^: 



(7.8) 



(7.9) 



(7.10) 



8=1 



which appear in the following expression for the microscopic unitary kernel (from eq. (|5.2j )): 

p(iV /+ l,a) (Ci;/X) ^) 



^ (2) (Ci,C 2 



a N f +i IM) 



CfC 2 Q - 



(-l)^Ai o) (Ci;M,>/=5 



A«(/x)(Ci-C2) v / nft(Ci + ^)(c 2 + /i 



(7.11) 
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The microscopic limit of the ipi functions is: 



Uo = f^ t i = 1 ,... tJV/ + 1 



AW( M , ( 7 i )(c+ ( 7?) v /n;L / 1 (c+^; 
AH(^,)v / n;= / i(c + ^ 2 ) v 



4>2N f +2{() 



AM(C,/*) 



A(- i )(M)vnSc+^) 



-a-2 



where <7j = for z = 1, . . . , Nf and <JN f +i = 0. The microscopic limit of the diagonal 
elements of the matrix A is (from eq. ( |6.4| )): 



2A t = - X {(3)dT 1 > a \»)P^ a \-vhv>), ^ = l,...,N f , (7.12) 
2A Nf+1 = -ax(/3)4? + + i 1 ' a " 1 V ) 0)P ( ^' a) (0;M) 

= " a * (/?) A%j A(-)(0 )M )- (7 ' 13) 

In equation ( [7.1 2|) the orthogonal polynomials P^ N f^ are evaluated at £ = —/if. As 
we have already discussed for the degenerate massive case, this situation should be con- 
sidered in a limiting sense of eq. ( 77?]) , that is C — — yu?. In particular, in the ratio 



A^(C, A*)/rii=i(C + one substitutes each element of the first line of the matrix in 
eq. ( ffT8|) with its derivative with respect to (, evaluated at ( = —/if. 
Finally, the elements of the matrix B and the terms iipi are given in eq. Q3.20D and 
eq. (|3.22|) , respectively. Let By and eip i be their microscopic large- N limit, respectively. 
By means of eq. ( |3.21| ), one obtains also C, Cqo and Cq. Putting together all of the terms 
above, we obtain the microscopic limit of the scalar kernel S^\x,y). It is very important 

(2) 

at this point to note that the number of corrections to the unitary kernel K N (x,y) ap- 
pearing in sffl (x, y) is indeed independent of N, because it depends only on the potential 
V(x), and thus in the double-microscopic large- iV limit the scalar kernel Sj^\x,y) is still 
obtained through the formulas ( |3.23|) , ( |3.24j ). Namely, in the orthogonal case (3 — 1, the 
scalar kernel is: 

2n 

-S (1) (Ci,C 2 ) = ^ (2) (Ci,C2)- £ [ACil-BACy^MC^e^) . (7.14) 

i<n ,j=l 
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In fact, the required substitution c — > 2c, which effectively means (3 — > 2/3 in this case, is 
compensated by the /3-scaling required for recovering the general case 1, as discussed at 
the beginning of this section. Moreover we have to put a — > 2a and Nf — > 2Nf, the latter 
meaning that all the masses are double degenerate. Indeed the degenerate-masses case 
does not present any peculiar problem, and by applying the same technique as described 



after eq. ( |4.9|) , the only consequence is that one has to consider higher-order derivatives of 
Bessel functions in all the determinants. 

In the symplectic case (3 = 4, one has to apply the substitution N — > 2N. Accordingly, 
the microscopic limit of the unitary kernel reads: 



lim 

N^oc 



N K, 



(2) 
2N 



•iL 

■2N'' 



2N' 



^ (2) (2Ci,2C 2 ), 



and the scalar kernel is 



£ (4) (Ci,C 2 



2d 



(2Ci,2G) -2 £ [AoCtfCok V><(2Ci) 6^(2$, 



i>n,j=l 



(7.15) 



(7.16) 



where once again the final /^-scaling is for recovering c = N(3 /2 ^ 1 . The only elements in 
the previous formulas which require further attention are the large-iV quantities B^ and 
eipj. In order to obtain them, we need to evaluate the large- N behaviour of the integrals 
defining the matrix elements Bij and eipi. When doing that, it could be advantageous to 
exchange the integrals with the large- iV microscopic limit, i.e. substituting the microscopic 
expansion of the functions ipi into the integrals. Such an interchange can be done under 



suitable smoothness assumptions. In |13| this issue has been investigated carefully for the 
massless case. It turns out that for /3 = 4 the procedure holds without any problems, 
whereas for (3 = 1 one has to consider also an additional contribution coming from the 
soft-edge of the spectrum. Since the case of one mass taken to zero or to infinity in our 
formulas reproduces results with one additional massless flavor and one massless flavor less, 
respectively, we shall assume that the interchanging procedure, when (3 = 4, is permitted 
in the intermediate mass region as well. Thus, in the (3 = 4 case we can write 



B^ 



OO /"OO 



JO 

oo 



dxdye(x - y) ^(x) ipj(y) 
dxe{x — y)ipj(x) , 



(7.17) 
(7.18) 



with i,j = 1, . . . , 2Nf + 2. The question whether such a simplification is allowed for (3 = 1 
case, can be addressed by means of numerical investigations. Nevertheless, with or without 
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exchanging the integrals and limits, from the very general result in eq. ( [7. 14 ) and eq. ( 7.16| ) 
one obtains analytic expressions for the microscopic limit of all the correlation functions 
of the general massive chiral ensembles for (3 = 1 and (3 = 4. For instance, the microscopic 
spectral density is the simplest one and it is given by p^(x) = S^'(x, x). Switching to the 
real eigenvalues £j of the Dirac operator, by means of £ 2 = x, one obtains the microscopic 
spectral density of the Dirac operator p^ (0 — 2|£|p^(£ 2 ).Q If universality arguments ap- 
ply to correlation functions of the general massive Chiral Ensembles, then the expressions 
obtained from eq. (|7.14|) and eq. (|7.16|) , through eq. (ft.ip, are the very natural candidates. 



8 Conclusions 

In this paper we considered the problem of computing the correlation functions for massive 
Dirac spectra of four-dimensional QCD. Starting from the fact that RMT is just a simple 
and effective tool for calculating actual spectral correlations in the infrared regime (which 
in principle also can be obtained from direct calculations in terms of finite-volume parti- 



tion function as in [34]] ) we consider a matrix model with matrices either in the %GOE 
or in the xGSE. The former corresponds to four-dimensional QCD with Nf fermions in 
the fundamental representation and 577(2) gauge theory, the latter corresponds to four- 
dimensional QCD with Nf fermions in the adjoint representation and SU(N C > 2) gauge 
theory. Matrix models with orthogonal and symplectic matrices naturally leads to the 
application of skew-orthogonal polynomials, which have the drawback that at present they 
are difficult to determine and handle in the microscopic large- N limit. Therefore we in- 
vestigated whether the recent technique proposed by H. Widom for dealing with OE and 
SE using standard orthogonal polynomials of the UE, could be applied effectively in our 
actual case. We have succeeded in deriving the orthonormal polynomials for the general 
massive fermion case, in an explicit and closed form. We found explicit formulas for the 
ij) functions which are the basic ingredients for computing the scalar kernels (x,y) for 
massive fermions, when (3 = 1 and (3 = 4. In particular, such scalar kernels are expressed 
as the unitary kernel Kjp (x, y) plus a finite number of corrections, which depend only 
on one-dimensional and two-dimensional integrals involving the functions fy. Cases with 
degenerate masses, or massless fermions can then be obtained just as simple limits of the 
formulas we have derived. We obtained fully analytical formulas for the scalar kernels 



9 In general, the fc-point correlation function flffi (£i, ...,£&) of the spectrum of the Dirac operator V is 
given by the change of variables f| = Xj in eq. (3.1). 
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and in the double microscopic large- N limit. From such formulas one can in princi- 
ple derive all the microscopic correlation functions of the Chiral Orthogonal and Chiral 
Symplectic Ensembles, with an arbitrary number of flavors, arbitrary masses and arbitrary 
topological charge. 

We emphasize that the issue of universality has not been considered in the present paper. 
However, the fact that the method of Widom can be applied to four-dimensional QCD 
with massive fermions is indeed encouraging and provides a new general framework where 
it might be possible to analyze some of the still open issues. We believe that the present 
framework is the most suitable setting for proving universality of these massive cases. 
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A Normalization factors 

The orthonormality condition of polynomials is an essential requirement for the proper 
use of Widom's method. Therefore, an explicit expression for the normalization factors of 
the orthogonal polynomials for the general Ay-flavor case is of primary importance. Even 
though it is a classical result, we choose to present in this appendix a derivation of the 
normalization factors h^ 1 '^, and of the highest coefficients kn f , i.e. eq. (|4.S| ) and eq. 
(fOl), respectively. Substituting Christoffel's formula eq. ( fOD into the orthonormality 
condition eq. (|4.6|) evaluated at i — j = n, yields: 

y/h!F f ' a) (m) = J + °° dxx a e- x Pi N f' a \x;m)det[A { n N f> a \x) 

= f °° dz^e-pW >.«)( x; m ) £ P£&\x)(-l) h det 

where we in the second equation used the so-called Laplace expansion of the determinant 
of a matrix with respect to its first row. Since / +o ° dx x a e~ x x n Pj°' a ^ (x) = 8 n j/k^ a \ for 
n < 3i one has that the only no n- vanishing term in the sum in eq. ( |A.1|) is when k = 0. 



A 



(N f ,a) 
n,k+l 



(A.1) 
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Therefore we obtain: 



^' Q) (m) 



Ill ) 



(0,o) 



det 



(N f ,a) 



1 Vi 



(A.2) 



with the coefficient /c^ ' ) given in eq. ( |4.5|) . The still unknown highest coefficients k^n 1 ^ 
can easily be determined as follows: 



fcCN>.«)( m ) = i im . 



,(JV/.«), 



E£oPS?(^)(-l)*det 



x' 



lim 



/in V/ ' a ''(m) + m. 



*£$,(-!)"' det 



A 



(N f ,a) 
n,N f +l 



(N f ,a) 



(A.3) 



111 ) 



In the first line we used Christoffel's formula eq. (177) and the Laplace expansion of the 
determinant, whereas the second line comes from the observation that only the last term 
in the sum contributes in the large- x limit. Combining eq. (|A.3| ) with eq. ( |A.2j ) one finally 
has 



e^(m) = (-l)^^||det 



A (*>.«)' 



det 



A 



(N f ,a) ■ 
n,N f +l 



111 ) 



\ 



(-ir> k}!'">k';i%- 



det 


a I 

yv n,iV / +l 


det 





where the coefficients k^' a ^ are given in eq. ( |4.5| ). 



(A.4) 
(A.5) 



B Some integrals 

In this appendix we evaluate the integrals in eq. (|5.9|) and eq. fl5.10|) , that is: 



T (N,j) 



\-fjV \Z, m J; 1 /N f ,a 



(Pi 



(N f ,a) 
N I 



z+m'j 



) Nf , a for % = 1,2, iV/ 



(B.1) 
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Using eq. (|4.16|) (i.e. adding m^f+i = 0) the integral I^'^ can also be written as 

D (JV>+l,a-l) / n N 

T (N,j) _ , p (JV)+i,a-i)/ n\ Li {z;m,0) 

J — VN {Z,m,()J, ) N f +l,a-l ) 

z 

therefore it can be obtained from l\ N in eq. (|B.1| ), by applying the substitutions Nf — > 
iV/ + 1, a -»• a - 1, {m} -> {m, 0} and i-*N f + l. 

The integral l\ N ^ can be evaluated as follows. First we apply eq. ( (4.19|) and then we use 
the orthonormality condition (|4.6|), that is 



N j 



r=0 s=0 



min(A r j) 

= </' Q) (m)4?' Q) (m) £ [^(-m?;^)] 2 . (B.3) 

r=0 

Using eq. ( |4 . 1 9| ) again, we finally end up with 

T (N,j) _ d N, l { m ) d j,i l m J p (Nf,a) , _2.__\ 
h ~ ,(N f ,a) J—, ^(NJ)K m ,' m ^ 

a min(N,j),i\ ) 

= ^j^P^i-ml m) . (B.4) 
Therefore the integral Jq ' is 

r(iV,i) _ ,(iV /+ l,a-l) , , (^+1,0,-1), , 
J - "max(AfJ),Af J +ll m ; U J- r min(Af,j) l U > m ' U J 

= ^Lax(Ar',j),Af^+l( m ' 0)-f > min / (ATj-)(0; m ) • (B-5) 

A final comment: one could wonder about the meaning of the "degenerate" expressions 
Pn f+1 ' a \0; m, 0) and Pn' a \— mf ; m). But as in the case of degenerate masses, they 
should be understood in the limit sense x — > or x — > — m?, respectively. Such a limit is 
easily evaluated from eq. ( |4.7|) by subtracting the z-th row from the first row of the matrix 
An ,a ' (eq. (|4.3|)), and dividing it by the term (x + mf) stemming from the product in 
the denominator of eq. fl4.7|) . Therefore, the limit x — > —mf of eq. fl4.7|) is equivalent to 
substituting all the polynomials Pj 0,a \x) in the first row of An 1 , with the corresponding 
derivatives py' a > \x) evaluated at x — —mf. 
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C The functions ipj(x) 



In this appendix we find a compact and useful expression for the functions ip. We will 
show the remarkable fact that every ip function is proportional to a polynomial Pn Nf ' a \ 
Let us first consider the functions in the Hilbert space TC. Substituting the coefficients Di 
eq. ( |5.13 ) into the linear combination eq. Q5.12 ) we obtain: 



,(Nf+l,a), 



N-l 



P 



(N f ,a) 



N-l 



2 \ (N f +l,a) , x 



^w( N f' a \x) 



(CI) 



where U{ = rrii for i = 1, . . . , Nf, <JN f +i = and in the second line we used the formula 



( }4.17|) . From eq. ( 4.18 ) we obtain: 

^(x) = 4" / ' a) (m)-^- 1 (X;m ' ai) 



k 



(N f +l,a). 



Nf + 1 



(C.2) 



N-l 



Also the ip functions in the orthogonal Hilbert space Ti 1 - can be written in a simplified 
form. In fact, let us first notice that by isolating the last term of the sum in eq. ( |4.19| ) 
and using again eq. ( 4.19| ) on the remaining sum, one has: 



P { N Nf ' a \x;m) 



that is 



d 



(N f ,a) 
N.i 



(m) 



P, 



(N f ,a) l 



N-l 



X\ 111 1 



P^\-mh^P N 



(N f -l,a) i 



x; m^i 



(C.3) 



P 



{N f ,a), 



j(JV>.«), 



NA 



m 



N 



(m) 



D (N f ,a) 
N-l 1 



x: m ) 



j^f^) ( ^ p (N f -l, a ) ( 2._ , p (N f -l, a ) 



N 



N 



x; m^i) 
(C.4) 

From eq. ( |5.15| ) and ( |5.17| ) (or ( |5.18| )), we exactly obtain the l.h.s. of last formula, therefore 

^N f +l+i(x) 



N 



,(JV>,«) ( 



m 



U (JV>-1,«) 



(m^) + 



mf) 



(C.5) 
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for i = 1, . . . , Nf and 



(C.6) 



N 



II I ) 



for i = iV/ + 1. 

Finally, let us remark that such compact forms eq. ( |C5| ) and eq. ( |U.6| ) are quite useful and 
effective in calculating the microscopic limit of ijj functions. 
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